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We present 2D Monte Carlo radiative transfer simulations of flattened prestellar cores. We 
argue the importance of observing prestellar cores near the peak of their emission spectra, 
Q ■ and we point out observable characteristic features on isophotal maps of asymmetric cores at 

?-H ' FIR wavelengths that are indicative of the cores' density and temperature structure. These 

^j^ , features are on scales 1/5 to 1/3 of the overall core size, and so high resolution observations are 

Cu ■ needed to observe them. Finally, we discuss the importance of determining the temperatures 

of prestellar cores with high accuracy. 
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1 Introduction 

Prestellar cores are condensations in molecular clouds that are on the verge of collapse or already 
collapsing. Their study is important for constraining the initial conditions for star formation 
(see Andre et al. 2000). Prestellar cores have typical sizes 2000 — 15000 AU, and typical masses 
0.05 — 10 Mq. They are cold with typical dust temperatures 7 — 20 K. Their density profiles 
tend to be flat in the centre and to fall off as r~^ (n = 2 — 4) in the envelope. Many authors have 
used Bonnor-Ebert (BE) spheres to represent prestellar cores (e.g. Alves et al. 2001). However, 
studies of the shapes of cores indicate that they are best described as triaxial ellipsoids, and they 
are close to being oblate spheroids (Goodwin et al. 2002). Additionally, numerical simulations 
of the collapse of turbulent molecular clouds indicate that prestellar cores are transient, non- 
spherical systems, which when projected onto the plane of sky appear to mimic BE spheres 
(Ballesteros-Paredes et al. 2004). Information about prestellar cores comes from molecular line 
observations (e.g. NH3, CO), but also from continuum observations, where cores are seen in 
emission (FIR, submm and mm) or absorption (NIR). In this paper, we focus on modelling the 
continuum emission of prestellar cores. 



2 Asymmetric models of prestellar cores 

We use 2D models to represent slightly asymmetric prestellar cores (see Stamatellos et al. 2004). 
Our goal is to capture generically the different features we might hope to detect on the isophotal 
maps of such cores. We use a Plummer-like density profile (Plummer 1915), that is modified to 
include azimuthally symmetric departures from spherical symmetry. This profile is simple and 
captures the basic observed characteristics of prestellar cores. Here, we present a model of a 
slightly flattened core (disk-like asymmetry; Fig. QJ left), with density profile 
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where no is the density at the centre of the core, and ro is the extent of the central region 
in which the density is approximately uniform. The parameter A determines the equatorial- 
to-polar optical depth ratio e , i.e. the maximum optical depth from the centre to the surface 
of the core {9 = 90°), divided by the minimum optical depth from the centre to the surface 
of the core {9 = 0° and 9 = 180°). The parameter p determines how rapidly the optical 
depth from the centre to the surface rises with increasing 9, i.e. going from the north pole 
at 6^ = 0° to the equator at = 90°. We assume that the core has a spherical boundary at 
radius i?core = 2 x 10^ AU, and no = 10^ cm"^, ro = 2 x 10^ AU, ry = 2, p = 4, and e = 2.5. 
The simulations are performed using Phaethon, a 3D Monte Carlo radiative transfer code 
(Stamatellos & Whitworth 2003). 

3 Dust temperature profiles 

The dust temperature (Fig. d right) drops from around 17 K at the edge of the core to 7 K 
at the centre of the core (also see Zucconi et al. 2001; Evans et al. 2001; Stamatellos & 
Whitworth 2003). We find that the dust temperature inside a core with disk-like asymmetry 
is 9 dependent, similar to the results of Zucconi et al. (2001) and Goncalves et al. (2004). As 
expected, the equator of the core is colder than the poles. The difference in temperature between 
two points having the same distance r from the centre of the core but with different polar angles 
0, is larger for more asymmetric cores (see Stamatellos et al. 2004). This temperature difference 
will affect the appearance of the core at wavelengths shorter than or near the core peak emission. 
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Figure 1: Left: Density on the x = plane for a flattened asymmetric core with equatorial-to-polar optical 

depth ratio e = 2.5 and p = 4. We plot isopycnic contours every 10°"^cm~^. The central contour corresponds to 

n = 10^"^cm~^. Right: Temperature on the x = plane for the same model, calculated with a Monte Carlo RT 

simulation. We plot isothermal contours from 8 to 18 K, every 2 K. 



4 SEDs and isophotal maps 

The SED of the shghtly asymmetric core we examine, is the same at any viewing angle, because 
the core is opticahy thin to the radiation it emits (FIR and longer wavelengths). However, 
the isophotal maps do depend on the observer's viewing angle. Additionally, they depend on 
the wavelength of observation. Phaethon calculates images at any wavelength, and therefore 
provides a useful tool for direct comparison with observations, e.g. at MIR (ISOCAM), FIR 
(ISOPHOT) and submm/mm (SCUBA, IRAM) wavelengths. We focus on wavelengths near 
the peak of the core emission (150-250 /im; we choose 200 /im as a representative wavelength). 
At 200 fim the core appearance depends both on its temperature and its column density in the 
observer's direction. This interplay between core temperature and column density along the line 
of sight results in characteristic features on the images of the cores (see Fig. [2|. Such features 
include (i) the two intensity minima at almost symmetric positions relative to the centre of 
the core, on the image at 30°, and (ii) the two intensity maxima, again at symmetric positions 
relative to the centre of the core, on the image at 90°. 
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Figure 2: Isophotal maps at viewing angles 0°, 30° 60° and 90°, for a flattened prestellar core with e = 2.5 

and p = 4 at 200 /xm. We plot an isophotal contour at 5 MJy sr~^ and then from 60 to 110 MJy sr~^, every 

5 MJy sr~^. There are characteristic symmetric features due to core temperature and orientation with respect to 

the observer. (We note the axes (x, y) refer to the plane of sky as seen by the observer). 



We conclude that isophotal maps at 200 /xm contain detailed information, and sensitive, high 
resolution observations at 200 fim are helpful in constraining the core density and temperature 
structure and the orientation of the core with respect to the observer. In Fig. [Hi left, we present 
a perpendicular cut (at x = 0) of the core images shown in Fig. El We also plot the beam size 
of the ISOPHOT C-200 camera (90'', or 9000 AU for a core at 100 pc) and the beam size of the 
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Figure 3: Left: A perpendicular cut through the centre of the core images presented in Fig. E) (but also including 
the background radiation field). We also plot the beam size of ISOPHOT and the beam size of the upcoming 
Herschel PACS/SPIRE (assuming a core distance of 100 pc). Centre: The effect of the parent cloud on cores. 
Temperature profiles of a non-embedded core (dashed lines), and of a core at the centre of an ambient cloud with 
Ay = 4 (dotted lines), and Am = 13 (solid lines). The upper curve of each set of lines corresponds to the direction 
towards the pole of the core {0 = 0°), and the bottom curve to the direction towards the core equator {0 = 90°). 
The difference between the two curves is indicative of the temperature gradient. The core is colder when it is 
inside a thicker parent cloud, and the temperature gradient within the core is smaller. Right: Same as on the 
left, but for a core embedded in a molecular cloud with visual extinction A\ — 4. The characteristic features at 
different viewing angles are weaker than in the case of a non-embedded core, but still observable. 

upcoming (2007) Herschel (13^^ or 1300 AU for the 170 /xm band of PACS; 17^' or 1700 AU for 
the 250 iiva band of SPIRE). ISOPHOT's resolution is probably not good enough to detect the 
features mentioned above. Indeed, a search in the Kirk (2002) sample of ISOPHOT observations 
(also see Ward-Thompson et al. 2002) does not reveal any cores with such distinctive features. 
However, Herschel should, in principle, be able to detect such features in the future. 



5 The effect of the parent cloud 

Cores are generally embedded in molecular clouds, with visual optical depths ranging from 2-10 
(e.g. in Taurus) up to 40 (e.g. in p Oph). Due to the ambient cloud the radiation incident 
on a core embedded in the cloud is reduced in the UV and optical, and enhanced in the FIR 
and submm (Mathis et al. 1983). Previous radiative transfer calculations of spherical cores 
embedded at the centre of an ambient cloud (Stamatellos & Whitworth 2003), have shown that 
embedded cores are colder (T < 12 K) and that the temperature gradients inside these cores 
are smaller than in non-embedded cores. 

Here, we examine the same core as before when it is embedded in a uniform density ambient 
cloud with different visual extinctions Ay. The ambient cloud is illuminated by the interstellar 
radiation field. Relative to the non-embedded core, the core embedded in an ambient cloud 
with Ay = 4 is colder and has lower temperature gradient (Fig. [HI centre). The isophotal 
maps are similar to those of the non-embedded core, but the characteristic features are less 
pronounced. This is because the temperature gradient inside the core is smaller when the core is 
embedded. In Fig. El right, we present a perpendicular cut through the centre of the embedded 
core image. It is evident that the features are quite weak, but they have the same size as in 
the non-embedded core (Fig. El left), and they should be detectable with Herschel^ given an 
estimated rms sensitivity better than ^1 — 3 MJy sr~^ at 170-250 /xm. 

Thus, continuum observations near the peak of the core emission can be used to obtain 
information about the core density and temperature structure and orientation, even if the core 
is very embedded {Ay ^ 10). 



6 Application: Modelling L1544 

We model L1544 (core distance D = 140 pc) using a slightly flattened density profile (e = 2 and 
^ = 4). Our goal is to fit the SCUBA 850 /xm image (Fig. El left), the SED data points (Fig. H 
right), and to calculate the temperature profile in the core. To fit the observational data we 
choose the asymmetry parameter e, based on the 850 fivn image of the core. Then we vary the 
mass of the core (by adjusting the central core density and the core flattening radius), so as to 
fit the submm and mm SED data. Finally, we vary the extinction through the ambient cloud in 
order to fit the FIR SED data. 
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Figure 4: Left: 850 /xm SCUBA image of L1544 (from Kirk 2002). Right: SED of the core. The hne corresponds 
to the model, and the points to the observed SED (data taken from Kirk 2002). 
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Figure 5: Left: 850 /xm isophotal map of the model. The core is viewed at an angle = 60°. Right: Dust 
temperature profile of the core and the immediate ambient cloud. Different lines correspond to different directions 

(from ^ = 0°, upper line, to ^ = 90°, bottom line). 



The best-fit model is obtained using central density no = 3.4 x lO^cm"^, inner flattening 
radius ro = 0.015 pc, core extent i?core — 0.075 pc, and an ambient cloud of visual extinction 
^cioud _ 23.7. These values are in good agreement with observations of L1544 (e.g. Ward- 
Thompson et al. 2002, Kirk et al. 2004). The calculated luminosity emitted from the core is 
^core = 0.10 Lq, in agreement with observations (Kirk 2002). The model fits the SED data 
very well (Fig. ^ right). The model also reproduces reasonably well the central region of the 
L1544 850 /im map (cf. Figs ^ and [HI left). However, it does not reproduce the asymmetries in 



the outer parts of the core, which are signatures of triaxiahty. 

The temperature at the edge of the core is Tedge = 8.7 K and at the centre Tcentre — 7.5 K 
(Fig. [3 left). The dust temperature we calculate with the model is lower by 2 — 3 K than the 
temperature estimated by Kirk (2002) using FIR (90, 170 and 200 /xm) ISOPHOT observations 
{Tiso = 10.2toi K). However, ISO observations have difficulty distinguishing the core from the 
ambient cloud, and thus this larger temperature may be due to the presence of the hotter ambient 
cloud in the observing beam. As a result of overestimating the dust temperature, the core mass 
calculated by Kirk (2002) is underestimated. The core mass is Mcore ^ ^a/ (^A(2dust)^A) 
(e.g. Andre et al. 1999), where Fx is the flux at a submm or mm wavelength, and kx is the 
assumed dust opacity per unit mass. If the core temperature is overestimated by 2 K (10 K 
instead of 8 K) then using the 850 fim flux the core mass is underestimated by a factor of 1.6 
(S850/xm(10K) = 1.6 X S85o/xm(8K)). The core mass of the model is M = 2.5 Mq which is 
similar to what Kirk et al. (2002) calculated (M-^^ ^^ = 2.5 ± 1). However, we note that Kirk 
uses a dust opacity KsbOfim — 0.01 cm^g~^, which is smaller by a factor of 2 than the dust opacity 
we use (/^850/im = 0.02 cm^g"-*^). Thus, to consistently compare the mass we need to divide the 
Kirk (2002) value by 2. 

We conclude that accurate estimates of core temperatures are important when calculating 
core masses from submm and mm observations, since e.g. overestimating temperatures by even 
just 1-2 K can lead to underestimating core masses by a factor of 2. 

7 Conclusions 

Far-infrared continuum maps of prestellar cores reflect both the column density and temperature 
fleld along the line of sight, and thus contain complementary information to the mm continuum 
maps that mainly trace column density. The radiative transfer models presented here show 
that the effect of the combined dust temperature and column density along the line of sight 
is to produce characteristic features in the FIR intensity maps. These features are useful for 
constraining the conditions in prestellar cores, and are expected to be present in cores with high 
enough temperature gradients (^2 K). 
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